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^ ■ Star formation is mainly determined by the observation of Ha radiation which 

■ is related to the presence of short lived massive stars. Disc galaxies show a 
strong cutoff in Ha radiation at a certain galactocentric distance which has led 
^ ! to the conclusion that star formation is suppressed in the outer regions of disc 
galaxies. This is seemingly in contradiction to recent UV observations^ that 
^ ' imply disc galaxies to have star formation beyond the Ha cutolf and that the 
star-formation-surface density is linearly related to the underlying gas surface 
. density being shallower than derived from Ha luminosities^. In a galaxy- wide 

O ! formulation the clustered nature of star formation has recently led to the insight 

I that the total galactic Ha luminosity is non-linearly related to the galaxy-wide 
Ph| star formation rate^. Here we show that a local formulation of the concept of 
O ■ clustered star formation naturally leads to a steeper radial decrease of the Ha 
^ • surface luminosity than the star-formation-rate surface density in quantitative 

ci '. agreement with the observations, and that the observed Ha cutoff arises natu- 



rally. 



The integrated galactic initial mass function (IGIMF) describes the mass spectrum of 



00 , 

0> . all newly formed stars in a galaxy. The IGIMF is calculated by adding up all stars of all 
^ I newly formed star clusters^' ^ and falls off more steeply with increasing stellar masses for 
massive stars^ than the canonical initial mass function (IMF) in each star cluster due to the 
Q ■ combination of two effects: the masses of the young star clusters are distributed according to 
0> . the embedded cluster mass function (ECMF), for which the upper mass limit is a function 
^ ' of the total star formation rate^ (SFR), and the stellar upper mass limit of the IMF is a 
function of the total star cluster mass^. Consequently, the total fraction of massive stars 
^ ' and therefore the total Ha luminosity drops faster with decreasing SFR than the SFR'^. The 
■ IGIMF theory has already been shown to naturally lead to the observed mass-met allicity 
relation of galaxies® and has received recent empirical verification in a study of IMF variations 
among galaxies^. 

In order to construct a quantitative local IGIMF-theory we introduce the local embed- 
ded cluster mass function (LECMF), 

aECMp(Me.,X,,) = ^^£^, (1) 

which defines the number of newly formed star clusters with mass Mgd per unit area at the 
location x, y in a disc galaxy. Observations^^ of Galactic star forming regions show that 
this function is a single part power law, ^lecmf oc M~i , with an index of /5 = 2. The 
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smallest cluster mass^, Med^min = 5 M©, should form at any place in the galaxy, whereas the 
most massive star cluster, Mcci,max,ioc(2^i y)) which can form locally is expected to depend on 
the local gas density, i.e. how much material is locally available for star cluster formation. 
Observations^'^ show that the most massive star cluster, Meci,max, of the whole galaxy is a 
function of the total galactic star formation rate. To express the upper hmit of the LECMF 
in dependence on the local gas surface density we write 



where Tigi^{x,y) and Sgas,o are the gas density at the location x,y and at the origin, respec- 
tively. The local mass of all star clusters between the two mass limits is determined by the 
local star formation rate which is described by the Kennicutt-Schmidt-law^' T^sFR{x,y) = 
A T,^^{x,y) with N — 1.4 being the widely accepted value^, whereas N — 0.99 follows 
from recent UV observations^. As UV emission is a star formation tracer that is much less 
sensitive to the presence of OB stars than Ha emission, the true exponent must be much 
closer to the value derived from UV observations. Thus, we chose = 1. Interestingly, the 
high gas-density part of the Kennicutt-Schmidt-plot^ based on FIR observations also has a 
flatter slope oi N — 1.08 and galaxy evolution models suggest N not to exceed unity in order 
to reproduce observed radial density proflles of disc galaxies^^, confirming our choice. The 
mass spectrum of all newly formed stars per unit area (the local integrated galactic initial 
mass function, LIGIMF) is calculated by adding up all newly formed stars of all young star 
clusters, and the Ha surface density follows by adding up the Ha flux contributions of all 
newly formed stars. The newly formed stars in each young star cluster are distributed ac- 
cording to the invariant canonical initial mass function (IMF)^^'-*-^ with a fixed lower mass 
limit but an upper mass limit depending on the total cluster mass''. Young star chistcrs 
above ~3000 Mq have constant Ha-light-to-mass ratios whereas smaller clusters are increas- 
ingly Ha under-luminous^. With decreasing star formation rate surface density the upper 
mass limit of the LECMF decreases, and consequently the fraction of under-luminous star 
clusters increases. Thus UV and Ha scale differently with the star formation surface density, 
gas surface density or galacotentric radius, respectively. A detailed explanation how the Ha 
surface luminosity is calculated is given in the Supplementary Discussion. 

The LIGIMF-theory is next applied to a sample of disc galaxies^^ with meassured gas 
surface densities and Ha surface luminosities of H ll regions averaged over annuli at different 
galactocentric radii. It is known that ionising photons emitted by massive stars can escape 
from well defined H II regions and lead to recombinations and thus Ha radiation in the 
surrounding diffuse ionised gas (DIG)^^. Using Ha emission as a star formation tracer this 
kind of photon leakage has to be taken into account to get an estimate of the true star 
formation rate. The study^^ of NGC 247 and NGC 7793 allows to construct a correction 
procedure in order to obtain the total Ha surface luminosities from the surface luminosities 
of H II regions only (see Supplementary Discussion) . 

For a linear star-formation law (A^ = 1) as derived from UV observations^ the LIGIMF- 
theory predicts an Ha surface luminosity as a function of the gas surface density which is 





in full agreement with the observations (Fig. 1). Additionally, the radial Ha profile derived 
in the LIGIMF-theory matches the observations perfectly (Fig. 2). The concept of clustered 
star formation resolves the discrepancies between Ha and UV observations completely. 

At first sight it might be objected that the LlGlMF-theory contradicts observations 
of the UV sources in the outer disc of galaxies: 5%-10% of all clusters in the outer discs of 
galaxies detected in UV have associated Ha emission^'''. The age estimates of the UV knots 
range up to 400 Myr. Clusters with Ha emission have ages <20 Myr as they are powered 
by short lived massive stars and therefore 5% of all observed UV knots are expected to have 
associated Ha emission in agreement with observations. The LIGIMF-theory predicts an 
overabundance of under-luminous star clusters beyond the Ha cutoff and a smaller number 
ratio of Ha to non-Ha-emitting UV knots is expected. But under-luminous does not mean 
no Ha emission. In the LIGIMF theory each young UV cluster is an Ha source, too, but UV 
and Ha luminosity scale differently with the cluster mass. Thus, this finding^'^ is entirely 
consistent with the LIGIMF-theory. The observed UV knots in the outer disc of M83 are 
at any time systematically smaller than their counterparts in the inner disc^^ in agreement 
with the fundamental basics of the LIGIMF-theory. Even one outstanding massive young 
star cluster in the outer region of M83 does not contradict this principle but is expected 
from a statistical point of view (see Supplementary Discussion). Furthermore, the M83 
FUV luminosity function of outer-disc stellar complexes is steeper than for the inner-disc 
population^^. A similar trend is reported in NGC 628 for the Ha luminosity function of H ll 
regions^^. In the LIGIMF-theory inner disc LECMFs have higher upper mass limits than 
in the outer disc LECMFs. Integration of the LECMFs over the outer regions leads to a 
steeper ECMF of the outer disc than the resultant ECMF of the inner disc, indicating that 
outer disc star formation complexes are systematically smaller than the inner disc ones. This 
integration effect is fundamentally the same as the IGIMF being steeper for dwarf galaxies 
with low global star formation rates than for disc galaxies with high star formation rates^. 

Previously, the Ha cutoff has been explaincd^^ by a drop of the local gas density below 
a critical density determined by the stability condition of a thin isothermal disc^*^' where no 
star formation can occur. In contradiction to this explanation recent UV observations^ reveal 
star formation outside the Ha cutoff, and dwarf galaxies^^ show star formation although 
their average gas density is lower than the critical density. Indeed, it has been shown that 
in regions with densities lower than the critical value, star formation can be driven by other 
than thermal instabilities^^. It has been argued that H II regions powered by the same 
massive stars are larger in a thin environment, i.e. at large galactocentric radii, than in a 
dense one and thus identical H ll regions become fainter in the outer galaxy. Therefore, it has 
been concluded^^ that the Ha surface luminosity should drop faster than the star formation 
rate surface density. Indeed, the surface brightness of individual H II regions should be 
fainter in the outer galaxy. But the Ha surface density considered in star formation laws 
refers to the total Ha luminosity per unit area of the galaxy and not to the cross section 
of the H II region. Identically powered H ll regions contribute equally to the Ha surface 
luminosity independently of their location in a thin or a dense gas environment. Thus, 
the proposed solution^^ does neither explain the Ha cutoff nor the different slopes of the 



UV-based and Ha-based star formation laws. It has been shown recently that a required 
minimum column density for massive star formation might exist^^ implying star formation 
with no massive stars in low density environments. However, this model predicts a top-heavy 
IMF for cloud column densities much larger than this threshold, for which no observational 
evidence exists^^, and allows no quantitative linkage of Ha luminosity and the star formation 
rate. Contrary to this previously existing work, the LIGIMF theory is in excellent agreement 
with the observed radial Ha- and UV-luminosity profiles (Fig. 2) and the Kennicutt- Schmidt 
star formation law (Fig. 1), and also allows the determination of star formation rates even 
in Ha-faint galaxy regions. 
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Figure 1 | Ha-luminosity surface density versus total gas surface density. The Ha- 

luminosity surface density versus the total gas surface density observed for seven disc 
galaxies^^ averaged over annuh at different galactocentric radii is plotted (black squares) 
after correcting for photon leakage from H ll regions (see Supplementary Discussion). These 
galaxies have a mean star formation rate of SFR=6.9 Mq yr~^ (3.2 - 16.4 Mq yr"^)^'-*^^, a 
mean total gas mass of Mgas = 2.1 ■ 10^° Mq (0.6 - 3.6 • 10^° M©)^'^^ and a mean scale length 
of Td = 4.4 kpc (3.9 - 5.2 kpc)^^"^^. These mean values define our model standard disc galaxy. 
For a choice of 7 = | the LIGIMF-theory predicts an Sna-Egas relation which matches the 
observations excellently (sohd fine) . Note that the underlying true star- formation density as 
derived from UV observations^ is directly proportional to the gas surface density {N — 1) 
and is shown after converting it into an Ha surface luminosity using the wrong linear Ken- 
nicutt Ha-SFR relation^' (dashed line) and shows the expected Sno-Sgas relation based 
on the classical picture which is in disagreement with the observations. Furthermore, the 
Ha-luminosity surface density in the high-luminosity part (Sna > lO'^^'^ erg s^^ pc^^) de- 
pends, for the correct LIGIMF-theory, on the gas surface density with a power of 1.4 (dotted 
line, extrapolated to low Ha surface luminosity) in agreement with the classical Kennicutt- 
Schmidt slope of = 1.4. The LIGIMF-theory puts the hitherto inconsistent Ha and UV 
observations in perfect agreement with each other. Both the steeper high-luminosity slope 
of = 1.4 and the Ha cutoff at low gas densities are two simultaneous outcomes of the 
LIGIMF-theory. The thick dotted line which coincides with the thick sohd curve shows a 
fitting function of the LIGIMF- model (see Supplementary Discussion) . 




Figure 2 | Star formation surface density versus scaled galactocentric radius. Radial 
distribution of the star formation surface density of nine disc galaxies based on UV^ (thin 
dashed-dotted lines) and Ra^^ (thin dotted lines) observations that rely on a wrong linear 
conversion^^ between the corresponding Ha-luminosity surface density and star-formation- 
rate surface density after correction for photon leakage (see Supplementary Discussion). 
The galactocentric radius is in units of the Ha threshold radius'^°. Over-plotted is the true 
underlying star-formation-rate surface density of our standard disc galaxy (thick black line) 
as defined in Fig. 1 and the model Ha surface luminosity (thick dotted black line) converted 
into a star- format ion surface density using the same linear convcrsion^^''^°. The LIGIMF- 
theory thus naturally accounts for the discrepant Esfr value at a particular radius. 



Supplementary Discussion 



This Supplementary Discussion provides more information on i) how the Ha surface 
luminosity is calculated, ii) the correction for photon leakage from H II regions, iii) the 
local applicability of the IGIMF theory, iv) a discussion on the sizes of typical star forming 
regions, v) a discussion of the presence of star clusters more massive than Meci,max,ioc and vi) 
the fitting formula for the Kennicutt-Schmidt law in the LIGIMF-model. 



1) Calculation of the Ha; surface luminosity 



In order to construct a local IGIMF theory we start with the IGIMF theory'^' devel- 
oped for the entire galaxy, and transfer all quantities into their corresponding local surface 
densities. 



Thus, we first define the local embedded cluster mass function (LECMF), 

eLECMp(Mee.,X,y) = ^^£^, (3) 

which is the number of newly formed star clusters with mass M^d per unit area at the 
location x, y in a disc galaxy. Observations^" of Galactic star forming regions show that the 
LECMF is a single part power law, ^lecmf mJ, with an index of /3 = 2. The local mass 
surface density of newly formed stars is 

,max,loc 

St Ssfr(x, y) = / aECMF(Mecl) Med dMed , (4) 

,min 

where Ssfr is the star-formation surface density, 6t ~ 10 Myr is the time span required to 
populate the cluster mass function completely^, and Meci,min = 5 Mq is the smallest cluster 
mass^. 



The observed most massive embedded star cluster in a galaxy is determined through 
the total star formation rate^'^, 

^^^^ = 84793 (^!:5_y '' . (5) 
Me \Meyr-^J 

To express the upper mass limit of the LECMF in dependence on the local gas surface density 
we write 



Sgas(^) y) 



7 



Mecl,max,loc(a^,2/) = M^a,ma^ ' • (6) 



^gas,0 



The maximal gas-mass surface density, associated with the position of the most massive 
embedded star cluster, in a disc galaxy, which has a total gas mass Mgas and a single- 
exponential gas disc, 

^gas{x,y) = Sgas,0 e'"'^'''' , (7) 



with a scale length r^, is determined by 



One may raise the objection if the crude assumption of an exponential gas disc model 
is sufficient to construct a standard galaxy as the outer galaxy discs may also be described 
by other expressions as for example a power law. The key issue of the theory developed here 
is that the gas surface density generally tends to decrease with increasing galactocentric 
radius. The main aim here is not the detailed modelling of individual galaxies but to ex- 
plain the characteristic discrepancies between Ha and UV observation and their traditional 
interpretations: i) a distinct star formation cutoff inferred from the Ha cutoff contrary to 
the extended radial star formation inferred from UV. ii) A steeper star-formation-law slope 
inferred from Ha (N=1.4) observations than from UV (N=1.0) observations within the Ha 
cutoff. The explanation through the present LIGIMF theory relies on the description of the 
gas surface density beyond the Ha cutoff. If the very outer gas disc is better described by 
for example a power law rather than by a single exponential law then this would not change 
the outcome of the LIGIMF theory. The radial Ha profiles of the disc galaxies examined 
in UV^ and Fig. 2 in the main Text are from ref. 30. Figure 9 of ref. 30 shows the total 
gas surface densities of some of these galaxies. As this is a log-lin diagram it can easily be 
seen that the general trend of how the gas surface density scales with radius can be well 
described by an exponential law up to ~2 Ha radii. If the aim is to construct detailed mass 
models of individual galaxies then an exponential law is inaccurate, but by far sufficient for 
our purpose here. 



The star-formation surface density is described by the Kennicutt-Schmidt-law^'^^, 

Esfr(x,?/)=AEJ1,(x,?/), (9) 

with N = 1.4 being the widely accepted value^, whereas N = 0.99 follows from UV 
observations'*^. As UV emission is a star formation tracer that is much less sensitive to 
the presence of OB stars than Ha emission, the true exponent N must be much closer to 
the value derived from UV observations. Thus, wc chose N = 1. Interestingly, the high 
gas-density part of the Kennicutt-Schmidt-plot^ based on FIR observations also has a flatter 
slope of = 1.08, and galaxy evolution models suggest A^ not to exceed unity in order to 
reproduce observed radial density profiles of disc galaxies^^, confirming our choice. 

The factor A is determined by integrating Supplementary equation (9) over the whole 

disc, 

SFRn^ , , 

A = - — 2 ■ (10) 

2 TT Sgas^o 



In each embedded star cluster the newly formed stars are distributed according to the 
invariant canonical initial mass function (IMF)^^'^"^ which is a two-part power law, m~"', with 



indices ai = 1.3 between miow = 0.1 and mi = 0.5 Mq , and 0:2 = 2.35 between mi = 0.5 Mq 

and mmaxi where mmax = ''^max(M2ci) is given by the maximum-stellar-mass-star-chister- 
mass relation^. This relation is well defined empirically and is a result of feedback-driven 
star formation on a star-cluster spatial scale of < few pc (Weidner, Kroupa & Goodwin, in 
prep.). The IMF in each star cluster, ^M^^p is normalised by the total star cluster mass, 

(M,,i) 

Meci = / m ^M,,i{m) dm . (11) 

Finally, the local integrated galactic initial mass function (LIGIMF) can be calculated by 
locally adding up all stars in all newly formed star clusters, 

CLiGiMF(m,a;,|/) = / Cm,,A^) Clecmf (M^chX^y) dM^a ■ (12) 

,min 

The surface density of ionising photons emitted by all new stars at the position x, y is 

/■mmax 

Nion,st{x,y) ^ CLiGiMF{rn,x,y) Nion,st{m) dm , (13) 

where A^ion,<5t(^) is the total number of ionising photons^ emitted by a star with mass m in 
time St. The Ha-luminosity surface density then follows from 

^Ha{x, y) = 3.02 X 10-'' erg Mon,5t(a:, y)/5t . (14) 

ii) Photon leakage 

In order to correct for the leakage of ionising radiation from H ll regions we use a 
comparative study'^ of the Ha radiation of well defined H II regions and their embedding 
diffuse ionised gas (DIG). This study of the two galaxies NGC 247 and NGC 7793 shows 
hat the radial and azimuthal surface luminosity of the DIG is highly correlated with bright 
H II regions and that the required power to sustain the DIG can only be met by the ionising 
radiation from massive star formation. In faint Ha luminosity regions the luminosity contri- 
bution of the DIG can exceed the one by distinct H II regions. We use the published values 
of the Ha surface luminosities of H II regions only and H II regions plus DIG in different 
regions of these galaxies to construct a correction formula. Supplementary Fig. 1 shows the 
fraction of the Ha luminosity due to H ll regions of the total (H ll + DIG) Ha luminosity 
as a function of the Ha luminosity in H ll regions. Two interpolating functions are used to 
describe the resultant relation. With 

V = ^<^!.Jy^^) (18) 
\ ^mi + DIG / 

and 

a; = logioSHii (16) 



the high luminosity region of NGC 7793 and all data points of NGC 247 can be easily 
interpolated (sohd line) by 

y{x) = 0.22 {x - 32.3) . (17) 
NGC 7793 alone can be interpolated (dashed hne) by 

y{x) = 0.22 {x - 32.3) - exp(-(a; - 28.2)) . (18) 

The total surface luminosities are then given by 

logio (^Hii + dig) ^x- y{x) . (19) 

Luminosities brighter than lO'^^'^ erg s^^ pc~^ need no correction. The influence of photon 
leakage from H ll regions for the Kennicutt-Schmidt diagram of Fig. 1 is shown in Sup- 
plementary Fig. 2. The corrected values (small black squares) are obtained by applying 
Supplementary equation (19) on the original values (open circles). The steeper [N — 1.4) 
Ena-Sgas relation left of the Ha cutoff and the position of the cutoff still persist. The lines 
are as described in the main text. Correcting for photon leakage therefore does not affect 
our results in any way. 

iii) Local applicability of the IGIMF theory 

Recently it has been demonstrated that the concept of clustered star formation has a 
significant influence on the appearance of galaxies in Ha^. The galaxy-wide Ha luminosity 
scales non-linearly with the underlying total star- format ion rate (SFR). As a first conse- 
quence it has been shown that a sample of the Sculptor group of dwarf irregular galaxies 
have much higher SFRs than previously deduced and that their SFRs scales approximately 
hnearly with their gas mass. The sample consists only of 11 galaxies and the vahdity of the 
obtained star formation relation might be questionable. But a subsequent investigation of a 
sample of 205 local star forming galaxies confirms this finding and reveals a strictly linear 
relation between the SFR of a galaxy and its total gas mass over five orders of magnitude 
in mass (Pflamm-Altenburg and Kroupa, in preparation). Furthermore, the local UV-based 
star formation law integrated over the disc of a galaxy is in perfect agreement with global 
Ha based and IGIMF corrected star formation law. Therefore, dwarf galaxies can be de- 
scribed very well with the IGIMF-theory, ie. by the concept of clustered star formation. For 
a radius of 3 kpc of a dwarf galaxy the reference area is about 28.3 kpc^. In this context 
"local" regions in a disc galaxy are annuli. E.g. an annulus between 10 and 11 kpc has an 
area of 66 kpc^, much larger than a dwarf galaxy for which the IGIMF theory is seemingly 
applicable. This means local regions in the context of disc galaxies have test areas larger 
than dwarf galaxies which we have mentioned to be excellently described by the IGIMF 
theory, i.e. by the concept of clustered star formation. 



iv) Sizes of outer disc star forming regions 

The basis of the LIGIMF theory is that X-UV discs result from local "stellar IMF bi- 
ases" . The Ha-luminosity — star-cluster-mass relation becomes non-linear^ for cluster masses 
< 3000 Mq. The LIGIMF theory predicts that inner-disc star-forming complexes sample 
the LECMF up to higher upper star-cluster mass limits than outer-disc star-forming com- 
plexes. As UV is a star formation tracer much less sensitive to the LIGIMF-effect (i.e. the 
supression of OB-star numbers in low-SFR regions), the UV/Ha luminosity ratio in the in- 
ner disc is smaller than in the outer disc. Thus, a basic requirement of the LIGIMF theory 
is that star forming regions in the inner disc are typically larger or more massive than the 
outer disc ones. A colour-magnitude diagram of UV knots observed in M83 differentiates 
into outer and inner disc UV knots (Fig. 3 in ref. 18). UV knots younger than 100 Myr 
have FUV-NUV colours (Fig. 2 in ref. 17). It can clearly be seen that the typical young 
(< 20 Myr) UV outer-disc knots of M83 have masses < 1000 Mq and the main bulk lies at 
K, 500 Mq while the inner knots extend to much higher masses. Only these star forming 
clusters (<20 Myr old) are of interest as only they contribute to Ha emission. Only one 
outer-disc object stands out with 10^ Mq. The presence of such a massive star cluster in the 
outer region is not forbidden statistically (see next Supplementary Section). Young inner- 
disc objects populate the CMD densely up to 5000 Mq and have masses which are typically 
significantly larger than outer-disc UV knots. Outer-disc knots which are much larger than 
1000 Mq are older than 20 Myr and can not be considered in the context of the Ha cutoff. 
But at these somehwat later ages the inner-disc knots are also more massive. Summarising: 
at any time the typical star forming regions are less massive in the outer disc than in the 
inner disc. Additionally, the median FUV luminosity of the inner-disc population is 3 times 
greater than in the outer region^^. And outer-disc complexes appear systematically fainter 
than their inner-disc counterparts^^. Furthermore, the M83 FUV luminosity function of 
outer-disc stellar complexes is steeper than for the inner-disc population. A similar trend 
is reported^^ for the Ha luminosity function of H ll regions in NGC 628. This is in full 
agreement with the concept of a varying upper mass limit of the LECMF [equation (2)]. 
If a galaxy is divided into two parts, an inner and an outer part, then the LEGMF of the 
respective part has to be integrated over the whole part of the galaxy to get the ECMF of 
each part separately. The inner part consists of regions where the LECMF has high local 
upper cluster mass limits. The outer part consists of regions where the LECMF has only 
low or intermediate local upper cluster mass limits. Thus, solving the surface integral over 
the LECMF it follows that the ECMF in the outer part is steeper than in the inner part. 
This is the same integration effect as for galaxies with a low SFR having steeper IGIMFs 
than galaxies with a high SFR^. The same holds for the FUV- luminosity and H ll-region 
luminosity function as observed. These observations are an indirect evidence that the up- 
per mass limit of the LECMF decreases towards galactic regions with lower star formation 
surface densities. 



Also the finding that 5%-10% of the UV knots in the extended disc of NGC 628 have 
associated H« emission^'^ is in-Unc with the LIGIMF theory. This ratio of Ha emitting 
sources to non-emitting sources is entirely reasonable given the age estimates of ^20 Myr 
for H II regions and Ri400 Myr for GALEX knots, i.e. 20/400 = 5 %. As explained above, 
an HcK-under- luminous star cluster is not the same as an Ha-non-luminous cluster. Each 
young UV source (<20 Myr) can appear as an Ha source. The crucial point is that the UV 
and Ha luminosities of star clusters of mass <3000 Mq scale differently with the star cluster 
mass^. Therefore, these observations^'' are in agreement with the LIGIMf theory. 

v) Local maximum cluster mass 

The presence, for example in outer-disc regions of galaxies, of outstanding star cluster 
masses is not forbidden, but expected. To demonstrate this we take the cluster mass proba- 
bility distribution density, /(M) = k M~^, between the lower. Mi, and an upper limit, M^, 
normalised by 

1= / /(M)dM. (20) 

J Ml 

N clusters are drawn from this distribution function with P — 2, Mi — 5 Mq, and Mu = 
10^ Mq, the mass of the most massive star cluster of this set is stored. This experiment is 
repeated several times. We then get a distribution of the most massive star cluster of a set 
of N star clusters. In general, the distribution function of the i-th massive star cluster of 
a set of star clusters can be constructed as follows: the probability that the i-th massive 
star cluster lies in the mass range from M to M + dM is /(M). There are (^) = A^ possible 
reahsations to choose the i-th massive star cluster out of A^. The probability that i — 1 more- 
massive star clusters lie in the mass range from M to M^ is [j^- f(M') dM'y . There 

are (^.T/) possible realisations to choose i — 1 clusters out of the remaining N — 1 clusters. 
Finally, the probability that N — i less-massive star clusters lie in the mass range from Mi to 

M is (^JI^^ f{M') dM'^ . There are (j^I*) = 1 possible realisations to choose N — i clusters 
out of the remaining N — i. Thus, the total distribution of the i-th massive star cluster is 
given by 

p,,^(M) = N (^^^~^^ ( f(M) dM^ /(M) /(M) dM^ . (21) 

For the most- massive star cluster {i — 1) and a single part power law for the mass function 
the distribution simplifies to 

/ k \N-1 

P^,n{M) = Nk ( ^(M^-/^ - M.'-f")] M-f , for /5 7^ 1. (22) 

The distribution of the most-massive star cluster is plotted in Supplementary Fig. 3 for 
A^ = 10,100,1000,10000 (drawing 10^ times from /(M)). Shown are the distributions 
obtained from Monte-Carlo simulations (grey histograms) and the analytical distributions 
(black curves) calculated using Supplementary equation (22). If a test area in the inner re- 
gion of the galaxy with a high star formation rate is specified one would expect, for example. 



1000 clusters per galactic test area and the expected most-massive star cluster and thus the 
upper mass limit of the LECMF in a typical inner region would be about 10^ Mq. In a test 
area of equal size in the outer regions were the star formation rate is lower by a factor of ten 
the expected number of star clusters per test area is reduced by a factor of ten, too. The 
expected most-massive star cluster would have a mass of shghtly less than 1000 M©, while 
a 10^ Mq star cluster is unhkely but possible. Thus, from the statistical point of view, the 
typical star forming regions in the outer galaxies have to be of lower mass than in the inner 
regions with a much higher star formation rate. 

vi) Fitting formula 

The relation between the theoretical Ha surface luminosity and the gas surface density 
can be very well fitted (Fig. 1) by a polynomial of fifth order, 

y = Go + aix + a2x'^ + a^x^ + a^x'^ + a^x^ , (23) 

where 

y = logio (Sgas/^oPC"^) , X = logio (^Ha/erg s"^ pc"^) . (24) 

The coefficients are listed in Supplementary Table 1. The underlying theoretical star forma- 
tion law [A^ = 1 in Supplementary equation (9)] is 

EsFR = 3.27 • 10-^° yr-i Sg^s • (25) 

By combining these two equations the LIGIMF theory allows a convenient and accurate 
conversion of an observed Ha luminosity surface density into a true star-formation-rate 
surface density. 
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Supplementary Table 1 | Coefficients for the Sna-Sgas fit- Coefficients of the fitting function 
[Supplementary equation (23)] of the relation between the Ha surface luminosity and gas surface 
density (Fig. 1). 
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Supplementary Figure 1 | Estimating the photon leakage from Hii regions. Ratio of the 
Ha-surface luminosity of H II regions only and the total Ha-surface luminosity (H ll regions 
plus DIG) versus Ha-surface luminosity of H ll regions only, based on a comparative study^^ 
of the two disc galaxies NGC 7793 and NGC 247 in annuli at different galactocentric radii. 




Supplementary Figure 2 | Influence of photon leakage on the derived Ha; surface den- 
sity. Ha-surface luminosity density corrected for photon leakage (small black squares) and 
uncorrected (open circles) versus gas surface density of the same galactic annuli. 
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Supplementary Figure 3 | Most-massive star cluster distribution. Distribution of the most 
massive star cluster of a set of N clusters from Monte-Carlo simulations (grey histograms) 
and the analytical treatment (black lines). 



